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A new approach to the electroweak properties of two- 
particle composite systems is developed. The approach is 
based on the use of the instant form of relativistic Hamilto- 
nian dynamics. The main novel feature of this approach is 
the new method of construction of the matrix element of the 
electroweak current operator. The electroweak current matrix 
element satisfies the relativistic covariance conditions and in 
the case of the electromagnetic current also the conservation 
law automatically. The properties of the system as well as the 
approximations are formulated in terms of form factors. The 
approach makes it possible to formulate relativistic impulse 
approximation in such a way that the Lorentz-covariance of 
the current is ensured. In the electromagnetic case the current 
conservation law is ensured, too. The results of the calcula- 
tions are unambiguous: they do not depend on the choice of 
the coordinate frame and on the choice of " good" components 
of the current as it takes place in the standard form of light- 
front dynamics. Our approach gives good results for the pion 
electromagnetic form factor in the whole range of momentum 
transfers available for experiments at present time, as well as 
for lepton decay constant of pion. 



I. INTRODUCTION. 

All atoms, nuclei, and the main part of the so called 
elementary particles are composite systems. That is why 
the constructing of correct quantitative methods of cal- 
culation for composite-particle structure is an important 
line of investigations in particle physics. In nonrelativis- 
tic dynamics there exist different correct methods which 
use model or phcnomenological interaction potentials. 
However, in the case of high energy one needs to develop 
relativistic methods. It is worth to note that now the ex- 
periments on accelerators, in particular at JLab are per- 
formed with such an accuracy that the treatment of tra- 
ditionally " nonrelativistic" systems (e.g. the deuteron) 
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needs to take into account relativistic effects. Relativis- 
tic effects are important also in the treatment of compos- 
ite systems of light quarks. In this case the relativistic 
effects are significant even at low energy. However, the 
relativistic treatment of hadron composite systems is a 
rather complicated problem. To solve it one needs, in 
fact, to solve a many-particle relativistic problem with, 
in addition, not precisely known interaction. Let us note 
that the use of the methods of the field theory in this 
case encounters serious difficulties. For example, it is 
well known that the perturbative QCD can not be used 
in the case of quark bound states (see, e.g., pLE|)- 

In the present paper we will use the relativistic con- 
stituent model which describes the hadron properties 
at quark level in terms of degrees of freedom of con- 
stituent quarks. The constituent quarks are consid- 
ered as extended objects , the internal characteristics of 
which (MSR, anomalous magnetic moments, form fac- 
tors) are parameters of model. As relativistic variant of 
constituent model we choose the method of relativistic 
Hamiltonian dynamics (RHD) (see, e.g., and the 

references therein). 

Our aim is to construct a relativistic invariant ap- 
proach to electroweak structure of two-particle compos- 
ite systems. The main problem here is the construction 
of the current operators @-|ll[ . 

It seems to us that RHD is the method the most ad- 
equate for our purpose. The use of RHD enables one to 
separate the main degrees of freedom and so to construct 
convenient models. 

We use one of the forms of RHD, namely the version 
of the instant form (IF) . 

Our approach has a number of features that distinguish 
it from other forms of dynamics and other approaches in 
the frames of IF. 

• The electroweak current matrix element satisfies 
the relativistic covariance conditions and in the 
case of the electromagnetic current also the con- 
servation law automatically. 

• We propose a modified impulse approximation 
(MIA). It is constructed in relativistic invariant 
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way. This means that our MIA does not depend on 
the choice of the coordinate frame, and this con- 
trasts principally with the "frame-dependent" im- 
pulse approximation usually used in instant form 
(IF) of dynamics, [j] 

• Our approach provides with correct and natural 
nonrelativistic limit ("the correspondence princi- 
ple" is fulfilled). 

• For composite systems (including the spin 1 case) 
the approach guarantees the uniqueness of the so- 
lution for form factors and it does not use such con- 
cepts as "good" and "bad" current components. 

• The approach describes correctly the spin Wigner 
rotation and this fact makes it possible to obtain 
the correct (QCD) asymptotic. 

The RHD method as the relativistic theory of com- 
posite systems is based on the direct realization of the 
Poincare algebra on the set of dynamical observables on 
the Hilbert space (see, e.g., |}p|| and the references 
therein) . RHD theory of particles lie between local field 
theoretic models and nonrelativistic quantum mechanical 
models. 

Let us describe briefly the main statements of RHD. 

As it is known (see, e.g., 12|), the relativistic invari- 
ance of a theory means that there exists (on the Hilbert 
space of states) the unitary representation of the inho- 
mogeneous group SX(2,C), which is the universal cov- 
ering of the Poincare group. The relativistic invariancc 
means the validity of the Poincare algebra commutation 
relations for the generators of space-time translations 
P M and rotations M^ v . To construct the representation 
SL(2, C) means to obtain these generators in terms of 
dynamical variables of the system. In the case of free 
particle system this program is not difficult to be re- 
alized and the generators P^ , Aif" have clear phys- 
ical meaning: P° = H — is the total energy opera- 
tor, P — (F 1 ,P 2 ,P 3 ) — is the operator of the total 
3-momentum, J = (M 23 , M 31 , M 12 ) — the operator of 

the total angular momentum, N = (M 01 , M 02 , M 03 ) - 
the generators of Lorentz boosts. 



1 It is known that correct impulse approximation (IA) real- 
ization in the frame of traditional version of IF dynamics en- 
counters difficulties: the standard IA depends on the choice of 
the coordinate frame. We show below that IA can be formu- 
lated in an invariant way, the composite system form factors 
being defined by the one-particle currents alone. 



In the case of interacting systems the situation is dif- 
ferent and the construction of the generators in terms 
of dynamical variables encounters some difficulties. To 
make clear the sense of these difficulties let us compare 
RHD with the nonrelativistic quantum mechanics. In 
nonrelativistic case it is the Galilean group that is the 
group of invariance. In the case of nonrelativistic inter- 
acting system the interaction operator enters (in addi- 
tive way) only the generator of time translations (energy 
operator). The interaction operator satisfies the usual 
conditions of invariance under translations and rotations 
and of independence on the choice of inertial coordinate 
frame. Under these conditions the algebraic relations for 
the generators of the Galilei group written in terms of dy- 
namical variables remain to be valid after the interaction 
including. This is the meaning of the Galilean invariancc 
of the theory. So, in nonrelativistic theory of interacting 
particles it is only one generator of Galilei group — the 
Hamiltonian — that contains the interaction. Other gen- 
erators have the same form as in the case of free particles. 
This way of interaction including in the observables alge- 
bra is unique in nonrelativistic case and gives the unique 
nonrelativistic dynamics — the dynamics governed by 
the Schrodinger equation. 

The situation in the relativistic case is quite different. 
The structure of the Poincare algebra is of such kind 
that the interaction including in the total energy operator 
alone results in the breaking of algebraic structure, i.e. 
in the relativistic invariance breaking. This is the conse- 
quence of the fact that the Lorentz transformations mix 
space and time. To preserve the algebra it is necessary 
to include the interaction in the other generators, too. 
Now the generators of Poincare algebra can be divided 
into two classes: containing interaction (interacting gen- 
erators) - hamiltonians, and without interaction (non- 
interacting generators). The latter present the so called 
kinematical subgroup. This division, that is the interac- 
tion inclusion in the Poincare group, is not unique. The 
different ways of such a division preserving the Poincare 
algebra result in different types of relativistic dynamics 
(see, e.g., [§). 

The idea of this approach — RHD — is originated by 
Dirac. In (l3) he considered different ways of descrip- 
tion of the evolution of classical relativistic systems - 
different forms of dynamics. Dirac separated the concept 
of time as a coordinate and that of time as a parameter 
defining the system evolution. Consequently, he defined 
three main forms of dynamics with different evolution pa- 
rameters: point (PF), instant (IF) and light-front (FF) 
dynamics. Each of these forms has its relative advan- 
tages and disadvantages. The total number of possible 
dynamics is actually five || so that the unique non- 
relativistic Hamilton description is changed in relativis- 
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tic case for, generally, five possibilities. Each of these 
possible dynamics is connected with a three-dimensional 
hypersurface in the four dimensional space. The initial 
conditions are given on these hypersurfaces and the evo- 
lution of the hypersurfaces is described. The kinemati- 
cal subgroups are the invariance groups of the hypersur- 
faces. In particular, these hypersurfaces are: hyperboloid 
i^i^o 2 , t > (PF), hyperplane t = (IF), light- 
cone surface x° + x 3 = (FF). 

It is worth to notice that RHD and the field theory 
are quite different approaches. The establishment of the 
connection between RHD and field theory is a difficult 
and as yet unresolved problem. Contrary to field theory, 
RHD is dealing with finite number of degrees of freedom 
from the very beginning. This is certainly a kind of a 
model approach. The preserving of the Poincare alge- 
bra ensures the relativistic invariance (see for details the 
Sec. II). So, the covariance of the description in the frame 
of RHD is due to the existence of the unique unitary 
representation of the inhomogeneous group SL(2, C) on 
the Hilbert space of composite system states with finite 
number of degrees of freedom. The success of composite 
models shows the validity of approximate relativistic in- 
variant description with fixed number of particles and a 
finite number of degrees of freedom. The similar situa- 
tion (as was pointed out in ||) takes place in the solid 
state theory: many of solid state properties are connected 
directly with its symmetry group and the actual form of 
the particle interaction plays less important role. RHD 
is based on the simultaneous action of two fundamental 
principles: relativistic invariance and Hamiltonian prin- 
ciple — and presents the tool the most adequate to treat 
the systems with finite number of degrees of freedom. 
It is worth to notice that the mathematics of RHD is 
similar to that of nonrelativistic quantum mechanics and 
permits to assimilate the sophisticated methods of phc- 
nomenological potentials and can be generalized to de- 
scribe three or more particles. 

Now the problem of the choice of the actual form of 
RHD arises. 

Some time ago it was proved that S'-matrices are 
equivalent in the different dynamics forms fl4|| . This fact 
is interesting but it does not mean the absolute equiva- 
lence of the forms. First, there are problems which can 
not be reduced to S'-matrix, e.g., the calculation of form 
factors. Second, one has to keep in mind that any con- 
crete calculation uses some approximations; the approx- 
imations usually used in different forms of dynamics are 
nonequivalent. 

Our point of view is the following. One must choose 
the form of dynamics adequate to the problem in question 
and to the approximations to be done. It seems us that 
this is in the spirit of RHD - the choosing of the adequate 



degrees of freedom. 

RHD is widely used in the theory of electromag- 
netic properties of composite quark and nucleon systems 
§, g, @, ||-||. It was shown that RHD not only 
presents an interesting relativistic model but is a fruitful 
tool and can compete with other approaches in describing 
the existing experimental data ( especially at small and 
moderate momentum transfers). At present time the FF 
dynamics is the most developed and used for composite 
systems g,[l]Jl||l|,|l|,|l^l . In fact, the light front dynam- 
ics has obvious advantages: a) the minimal number of 
operators containing the interaction (three); b) the sim- 
ple relativistic invariant separation of the variables into 
classes of " internal" and " external" variables (while con- 
sidering the approach as Hamiltonian theory with fixed 
particle number); c) the simple vacuum structure for the 
light-front perturbative field theory. However, there are 
some difficulties in the FF RHD approach. In particular, 
it was shown |l5|,^9| that the calculated electromagnetic 
form factors for the systems with the total angular mo- 
mentum J = 1 (the deuteron, the p - meson) vary signif- 
icantly with the rotation of the coordinate frame. This 
ambiguity is caused by the breaking of the so called angle 
condition (l^j2^] , that is really by the breaking of the 
rotation invariance of the theory. Some of the difficulties 
of FF dynamics are discussed in [^0|. A possible way 
to solve the problem by adding some new (nonphysical) 
form factors to the electromagnetic current was proposed 
(see |3l| and references therein). 

A different approach to the problem was proposed re- 
cently in Ref. [ill , where a new method of construction 
of electromagnetic current operators in the frame of FF 
dynamics was given. The method of |Tl[| gives unam- 
biguous deuteron form factors. However, as the authors 
of [ [TT| note themselves, their current operator and the 
one used in Ref. || are different, since both of them 
are obtained from the free one, but in different reference 
frames, related by an interaction dependent rotation. 

All these facts naturally cause authors to consider 
other forms of relativistic Dirac dynamics. 

Recently the PF dynamics was considered in the pa- 
pers [| 23 2^^]. The authors used PF dynamics to 
calculate the processes under investigations in JLab ex- 
periments. It is worth to notice that these experiments 
enlarge the interest to the RHD approach - the relativis- 
tic theory which can be used in the region of soft pro- 
cesses. 

Now we present a relativistic treatment of the prob- 
lem of soft electroweak structure in the framework of 
another form - IF of RHD. IF of relativistic dynamics, 
although not widely used, has some advantages. The cal- 
culations can be performed in a natural straightforward 
way without special coordinates. IF is particularly con- 
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venient to discuss the nonrelativistic limit of relativistic 
results. This approach is obviously rotational invariant, 
so IF is the most suitable for spin problems. 

We describe the dynamics of composite systems (the 
constituent interaction) in the frame of general RHD ax- 
iomatics. However, our approach differs from the tra- 
ditional RHD by the way of constructing of matrix ele- 
ments of local operators. In particular, our method of 
description of the electromagnetic structure of compos- 
ite systems permits the construction of current matrix 
elements satisfying the Lorentz-covariance condition and 
the current conservation law. 

To construct the current operator in the frame of IF 
RHD we use the general method of relativistic invariant 
parameterization of matrix elements of local operators 
proposed as long ago as in 1963 by Cheshkov and Shi- 
rokov §§. 

The method of |3^] gives matrix elements of the op- 
erators of arbitrary tensor dimension (Lorentz-scalar, 
Lorentz- vector, Lorentz-tensor) in terms of a finite num- 
ber of relativistic invariant functions - form factors. The 
form factors contain all the dynamical information on the 
transitions defined by the operator. That is why a system 
can be described in terms of form factors. 

The method of parameterization is similar in spirit to 
the method of presentation of matrix elements of irre- 
ducible tensor operators on the rotation group in terms of 
reduced matrix elements. This method extracts from the 
matrix element of a tensor operator a part defining sym- 
metry properties and selection rules following the well 
known Wigner— Eckart theorem. 

In the review Gj two possible variants of such kind of 
representation of matrix elements in terms of form fac- 
tors are presented - the elementary-particle parameteri- 
zation and the multipole parameterization. The variant 
of parameterization given in ]3^ | is an alternative one. 
In p^] the authors propose the construction of matrix 
elements in canonical basis so it can be called canonical 
parameterization. This method was developed for the 
case of composite systems in jj3|,[34). The composite- 
system form factors in this approach are in general case 
the distributions (generalized functions), they are defined 
by continuous linear functionals on a space of test func- 
tions. Thus, for example, the current matrix elements 
for composite systems are functionals, generated by some 
Lorentz-covariant distributions, and the form factors are 
functionals generated by regular Lorentz-invariant gen- 
eralized functions. We demonstrate these facts below, in 
Sec. Ill, using a simple model as an example. 

It is worth to notice that the statement that the form 
factors of a composite system are generalized functions is 
not something exotic. This fact takes place in the stan- 
dard nonrelativistic potential theory, too (see Sec.III(F)). 



The use of canonical parameterization permits to de- 
scribe the electroweak properties of composite systems 
satisfying the Lorentz-covariance condition on each stage 
of calculation and to satisfy the electromagnetic cur- 
rent conservation law when describing the electromag- 
netic properties. In our formalism it is necessary to for- 
mulate the composite model features in slightly unusual 
way in terms of matrix elements which are generalized 
functions. 

In particular, the relativistic impulse approximation 
(IA) has to be reformulated. 

Let us remind the physics of IA. In IA a test particle 
interacts mainly with each component separately, that 
is the electromagnetic current of the composite system 
can be described in terms of one-particle currents. In 
fact, the composite-system current is approximated by 
the corresponding free-system current. This means that 
exchange currents are neglected, or, in other words, that 
there is no three-particle forces in the interaction of a 
test particle with constituents. It is well known that 
the traditional IA breaks the Lorentz-covariance of the 
composite-system current and the conservation law for 
the electromagnetic current (see, e.g., M for details). As 
we show in the Sec.III(C) one can overcome these diffi- 
culties if one formulates IA in terms of form factors. 

It is worth to notice that all known approaches (in- 
cluding the perturbative quantum field theory (QFT)) 
encounter difficulties while constructing a composite- 
system current operator satisfying Lorentz-covariance 
and conservation conditions. This problem is now dis- 
cussed widely in the literature |Q-[TTJ . To satisfy the con- 
servation law in the frame of Bethe-Salpeter equation 
and quasipotential equations, for example, it is necessary 
to go beyond IA: one has to add the so called two-particle 
currents to the current operator. In the case of nucleon 
composite systems these currents are interpreted as me- 
son exchange currents In the case of deuteron this 
means the simultaneous interaction of virtual 7- quanta 
with proton and neutron. However, in Ref. p5fl it is 
shown that the current conservation law can be satisfied 
without such processes, although they contribute to the 
deuteron form factor. It seems that at the present time 
there is an intention to formulate IA with transformed 
conservation properties without dynamical contribution 
of exchange currents |ll],|2(||3l| . 

Our formalism also gives, in fact, the description of 
the covariance properties of the operators in terms of 
many-particle as well as one-particle currents. However, 
the important feature of our formalism is the fact that 
form factors or reduced matrix elements describing the 
dynamics of transitions contain only the contributions of 
one-particle currents. 

So, our approach to the construction of the current 
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operator includes the following main points: 

1. We extract from the current matrix element of com- 
posite system the reduced matrix elements (form factors) 
containing the dynamical information on the process. In 
general these form factors are generalized functions. 

2. Along with form factors we extract from the matrix 
element a part which defines the symmetry properties of 
the current: the transformation properties under Lorentz 
transformation, discrete symmetries, conservation laws 
etc. 

3. The physical approximations which are used to calcu- 
late the current are formulated not in terms of operators 
but in terms of form factors. 

In this paper we present the main points of our ap- 
proach. To make it transparent we consider here only 
simple systems with zero total angular momenta, so that 
technical details do not mask the essence of the method. 
We demonstrate the effectiveness of the approach by cal- 
culating the pion electroweak properties. In this case the 
canonical parameterization is very simple and can be re- 
alized without difficulties. The case of more complicated 
systems needs using a rather sophisticated mathematics 
for canonical parameterization of local operator matrix 
elements and will be considered elsewhere. 

The paper is organized as follows. In Sect. II we re- 
mind briefly the basic statements of RHD, especially of IF 
RHD. Two bases in the state space of composite system 
are considered: the basis with individual spins and mo- 
menta and the basis with separated center-of-mass mo- 
tion. The Clebsh-Gordan decomposition for the Poincare 
group which connects these bases is given. The IF wave 
functions of composite systems are defined. In Sect. Ill 
our approach to relativistic theory of two-particle com- 
posite systems and their electroweak properties is pre- 
sented. A simple model is considered in details: two spin- 
less particles in the S-state of relative motion, one of the 
particles being uncharged. The electromagnetic form fac- 
tor of the system is derived. The standard conditions for 
the current operator are discussed. The modified impulse 
approximation (MIA) is proposed. The results of IA and 
MIA are compared. The nonrelativistic limit is consid- 
ered. The connection between the presented version of IF 
RHD and the dispersion relations in mass is established. 
In Sect. IV the developed formalism is used in the case 
of the system of two particles with spins 1/2. The pion 
electromagnetic form factor and the lepton decay con- 
stant are derived. The model parameters are discussed 
and the comparison of the results with the experimental 
data is given. The results of calculations in IA and MIA 
are compared and are shown to differ significantly. In 
Sect.V the conclusion is given. 



II. RELATIVISTIC HAMILTONIAN DYNAMICS. 

In this Section some basic equations of RHD and 
some relations from relativistic spin theory are briefly 
reviewed. 

The relativistic invariance of a theory means that the 
unitary representation of the Poincare group is realized 
on the Hilbert space of system states (see, e.g., p2|). 
In this case the structure of the Poincare algebra can be 
defined on the set of observables (M^ ,P a ) 

[M^,P a ] = -i(g^P v - fl^P"), 

[M^, M ap ] = -i(g^M vp - g va M» p ) - (a «-► p), 

[P' i ,P v ]=0. (1) 

In (Q) g^ is the metric tensor in Minkowski space. 
As we have mentioned, in the case of the particles with 
interaction the realization of the Poincare algebra on the 
set of observables is more complicated than in the case of 
the free particles. Let us consider one key commutator 
from the set of commutators ([!]) (see, e.g., Q): 

[p3N k ] = i 5 jk H (2) 

(we use notations given in the Introduction; j, k = 
1,2,3). 

Since H is interaction dependent for non-trivial sys- 
tems, either P , N, or some combination of P and N 
also must be interacting. To preserve the commutation 
relations (^) one has to make other generators depending 
on the interaction, too. So, the generators occur to fall 
into two groups: the generators which are independent of 
the interaction and form the so called kinematical sub- 
group, and the generators depending on the interaction 
- Hamiltonians. This division is not unique. Different 
ways to obtain kinematical subgroups result in different 
forms of dynamics. 

In this paper we use the so called instant form dynam- 
ics (IF). In this form the kinematical subgroup contains 
the generators of the group of rotations and translations 
in the three-dimensional Euclidean space: 

f, P. (3) 

The remaining generators are Hamiltonians (interaction 
depending): 

P° , N . (4) 

The additive including of interaction into the mass 
square operator (Bakamjian-Thomas procedure ]36[ ], see, 
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e.g., |jj for details) presents one of the possible technical 
ways to include interaction in the algebra (|l]): 



Ml + U 



(5) 



Here Mq is the operator of invariant mass for the free 
system and Mj - for the system with interaction. The 
interaction operator U has to satisfy the following com- 
mutation relations: 



P, U = 



J,U = 



Vp,U =0 



(6) 



These constraints ([}]) ensure that the algebraic relations 
(1) are fulfilled for interacting system. The constraints 
(6) are not too strong. For instance, a large class of non- 
relativistic potential satisfies @. The relations @ mean 
that the interaction potential does not depend on the to- 
tal momentum of the system. This fact is well established 
for a class of potential, for example, for separable poten- 
tials p^ . Nevertheless, the conditions (||) and (^) can 
be considered as the model ones. There exists another 
approach |38| where a potential depends on the total mo- 
mentum but that approach is out of scope of this paper. 

In RHD the wave function of the system of interacting 
particles is the eigenfunction of a complete set of com- 
muting operators. In IF this set is: 



M 



J 



Ja, P 



(7) 



j 2 is the operator of the square of the total angular mo- 
mentum. In IF the operators J 2 , J3 , P coincide with 
those for the free system. So, in (fj]) only the operator 
Mf depends on the interaction. 

To find the eigenfunctions for the system (Q) one has 
first to construct the adequate basis in the state space of 
composite system. In the case of two-particle system (for 
example, quark-antiquark system q q) the Hilbert space 
in RHD is the direct product of two one-particle Hilbert 
spaces: Ti q q = H q ® Ttq. 

As a basis in H qq - one can choose the following set of 
two-particle state vectors: 

I Pi , mi; P2 , m 2 ) = | pi mi) ® | pi m 2 ) , 

(p 1 m\p t m') = 2p S(p-p')5 mm > . (8) 

Here p\ , p 2 are 3-momenta of particles, mi , to 2 — spin 
projections on the axis z, po = yfp 1 + M 2 , M is the 
constituent mass. 

One can choose another basis where the motion of the 
two-particle center of mass is separated and where three 
operators of the set (Q) are diagonal: 

\P, >fs, J, I, S, mj), 



(P, A J, I, S, mj\P', V?, J', I', S\ mj,) 
N CG 5 {3) (P - P')S{^/~s - y/j)8jj.8w&ss'& mjmjl , 



N, 



CG 



(2Po) 
8k J 



, k = l^/s-AM 2 



(9) 



Here P^ — {p\ + pi)^, P?i — s, y/s is the invariant 
mass of the two-particle system, / — the orbital angu- 
lar momentum in the center-of-mass frame (C.M.S.), 
S 2 = (5i + S 2 ) 2 = S(S + 1) , S — the total spin in 
C.M.S., J — the total angular momentum with the pro- 
jection mj. 

The basis (||) is connected with the basis (^) through 
the Clebsh-Gordan (CG) decomposition for the Poincarc 
group. The decomposition of the direct product (||) 
of two irreducible representations of the Poincare group 
into irreducible representations ^ has the following form 



I Pi, mi; P2,m 2 , 



P, Vs, J, I, S, mj) 



x ( Jmj\S I m s mi )Y l * rni (1? , tp) (S m s |l/2 1/2 m x m 2 } 
x (mi|D 1 / 2 (P,pi)|mi)(m 2 | J D 1 / 2 (P,p 2 )|m 2 ). (10) 

Here the sum is over the variables mi, m 2 , to;, ms, I, 
S, J, mj. p = (pi — p 2 )/2, p — \p\, 1? , <p are the spheri- 
cal angles of the vector p in the C.M.S., Yi mi - a spher- 
ical harmonics (star means the complex conjugation), 
( S ms 1 1/2 1/2 mi to 2 ) and ( Jmj\S I ms mi } are the CG 
coefficients for the group SU(2), (m\ D 1 ^ 2 (P,p) \m) - the 
three-dimensional spin rotation matrix to be used for 
correct relativistic invariant spin addition. 

Let us discuss briefly the relativistic properties of spins. 
It is known that the Lorentz transformation for spins is 
momentum depending (see, e.g., Q). So, to perform 
Lorentz invariant spin addition for particles with differ- 
ent momenta p and p' one has to "shift" the spins to 
the frame where the momenta are equal to one another. 
The spin transforms following the so called small group 
which is isomorphic to rotation group and thus, the op- 
erator of such a " shift" is a 3-dimensional rotation matrix 
D(a,/3,j) . The Euler angles a, (3, 7 can be written in 
terms of the components of the vectors p and p' . In this 
way the "transplantation" of spins on one and the same 
momentum is realized. To understand what means this 
"transplantation" let us consider an example: one par- 
ticle has the momentum pi, mass M\ 1 spin j and spin 
projection to, while another particle with momentum p 2 
and mass M 2 has no spin. In the case of free particles 
the vector state of this system is 
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\px,M x ,j ,m;p 2 ,M 2 ) = \pi , M x ,j,m) ® \p 2 ,M 2 ). 

(11) 

The two-particle vector state in the case when it is the 
first particle that has no spin is 

| pi ,M X ; p 2 ,M 2 ,j,m) = \ pi ,M X ) ® \p 2 ,M 2 ,j ,m). 

(12) 

In nonrelativistic angular momentum theory the states 
( pi] ) and (12) are identical. They describe the two- 
particle system with momenta p\ and p 2 and total spin 
j with the projection m. In both cases the total spin 
can be obtained in simple way and is equal to the spin j 
of the first particle for (11) or of the second particle for 

In relativistic theory the states (O) and (O) differ es- 
sentially. The difference is caused by the fact that the 
states transform from one inertial coordinate system to 
another in different ways. As was mentioned before the 
Lorentz transformation for spin depends on the particle 
momentum and the spins in ( |ll| ) and ( |l2] ) correspond 
to particles with different momenta. In relativistic case 
the states ([ll]) and ( |l2| ) coincide only if the particles mo- 
menta are equal. In general case, to connect the state 
vectors ( pi] ) and ( p2| ) one has to "shift" the spin, for 
example, to shift the spin in (^) into the frame where 
the second particle has the momentum p\. This shift- 
ing transformation is realized by the matrix D (p 2 ,pi) 
which belongs to the small group. Let us consider the 
state vector 



| pi , Mi ) ® D 3 (p 2 , Pl ) | p 2 , M 2 , j ,m). 



(13) 



Let us remind that a transformation belonging to the 
small group does not act on momenta. It is easy to 
see that the resulting vector describes the same state as 
( pi] ) and transforms from one frame to another in the 
same way as ([ll]). To show this let us use the equation 
D(p' ,p") D{p" ,p) — D(p',p). Now the following covari- 
ant equality is valid: 

I Pi ,Mi ,j ,m; p 2 ,M 2 ) 
= J2 |pi,Mi;ft,M a ,j,m / )<m / |I> J '(p2,pi)|m). (14) 

m' 

As one can see from (|l4| ) the spin has "changed" the 
momentum. As this operation is very important for the 
understanding of the parameterization below let us for- 
mulate it in other words, namely, in terms of the gener- 
ators of the Lorentz transformations. 

The Lorentz-transformation generator for the state 
(O) is of the form: 



N =i P°^--^TM' ^ = ipo^- (15) 
op po + M op 

If we perform the transformation (03), the D- matrix 
transforms the generators in the following way: 

D j (p 2 , Pl ) (N, + N 2 ) [D*(p2,pi)]- 1 =N X + N 2 . (16) 

It is just the equations (|ll|) and ([!(]) (and only these 
equations) that present the exact statement that D- 
function shifts spin from one momentum to another one. 
Similar equations (however, more complicated) are valid 
in the case of two non-zero spins. In the case of spin 1/2 
the D - function fl34| has the form 



D 1/2 { Pl , P2 ) = cos (|)-2i(fcj) sin ( 



[PlP2_ 

[pipy 



lj = 2 arctan 



|[pip 2 ]| 



(piO + Mi)(P20 + M 2) - (PlP 2 ) 



(17) 



We will use it below. 

Let us make a remark concerning the invariance of the 
decomposition jTo|). The total spin S and the total or- 
bital angular momentum I in (|10| ) play the role of invari- 
ant parameters of degeneracy. However, the square of the 
total spin S = (Si + S 2 ) is not invariant. But one can 
define the total spin square in invariant way as follows: 

{[D^ipuP^S^D^ipuP)} 

+ [D^{p 2 ,P)}- 1 S 2 [D^(p 2 ,P)]y =S(S + l). (18) 

Here P is the center-of-mass momentum. One can see 
that in C.M.S. the definition ( p"8| ) coincides with the def- 
inition §). Similarly, one can define the orbital moment 
I in invariant way. 

The described spin rotation effect in ( ^0| ) is a purely 
relativistic effect. If one takes it into account, one obtains 
interesting observable effects [ B4[ . 

To obtain the basis vectors (g) in terms of vectors (||) 
one has to inverse (|l0|). The final equation has the form: 

IP, y/a, J, I, S, mj) 



Ei dpi dp 2 
I 7T— 7T— \Pi,mi, p 2 ,m 2i 



2pio 2p 
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x (pi ,mi; p 2 ,m 2 \P, y/s, J, I, S, mj) . (19) 



Here 



N = N 1 +N 2 



{p u m 1 ;p 2 ,m 2 \P, y/s, J, I, S, mj\ 
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^[X(s, M 2 , M 2 )]- 1 ' 2 2P 5(P - Pl - pa) 



{ mi| ^ 1/2 (pi P) |mi )< m 2 1 D x / 2 (p 2 P)\m 2 ] 



{l/2l/2m 1 fh 2 \Sm s )Y lmi {d,ip) 



J u f( k ) k 2 dk = l 



(23) 



Let us note that for composite c 
sometimes instead of equation (t 



uark systems one uses 
the following one: 



J uf(k)k 2 dk = 1 . 



(24) 



x (S lm s mi | Jmj) . 

Here X(a, b, c) — a 2 + b 2 + c 2 — 2(ab + bc + ac). To obtain 
( |l9| ) the decomposition in terms of spherical harmonics 
and the summation of all of the momenta to give the total 
momentum J were performed in the C.M.S. and then the 
obtained result was shifted to arbitrary coordinate frame 
by use of D- functions. 

It is on the vectors (S), ([l9]) that the Poincare-group 
representation is realized in the vector state space of two 
free particles. The vector in representation is determined 
by the eigenvalues of the complete commuting set of op- 
erators: 



M 2 



P 2 , J\ h 



(20) 



The parameters S and I (as was mentioned) play the role 
of invariant parameters of degeneracy. 

As in the basis (||) the operators J 2 , J3 , P in (Q) are 
diagonal, one needs to diagonalize only the operator Mf 
in order to obtain the system wave function. 

The eigenvalue problem for the operator M 2 in the ba- 
sis (^) has the form of nonrelativistic Schrodinger equa- 
tion (see, e.g., {§). 

The corresponding composite-particle wave function 
has the form 

(P', s/7, J', l', S',m'j\p c ) = 
= N C S(P' - fJSjj'Snjm'j <p$ s ,{tf) , (21) 



N c = V2pc 



N, 



CG 



4fc' 



Pc 



is an eigenvector of the set (|7|); J(J + 1) and mj are 
the eigenvalues of J 2 , J3, respectively (Eqs. (Q), (|2(j|)). 

The two-particle wave function of relative motion for 
equal masses and total angular momentum and total spin 
fixed is: 

i pf s (k(s)) = ^ui(k)k, (22) 
and the normalization condition has the form: 



Here n c - is the number of colours. The wave func- 
tion ( |22] ) coincides with that obtained by "minimal rela- 
tivization" in |39|. The normalization factors in ( |22] ) in 
this case correspond to the relativization obtained by the 
transformation to relativistic density of states 



k 2 dk 



k 2 dk 



(25) 



It is worth to notice that wave functions in RHD (for 



(Zlh defined as the eigen- 
) in general are not the 



example, the wave function 
function of the operators set ( 1 
same as relativistic covariant wave functions defined as 
solutions of wave equations or as the matrix elements of 
local Heisenberg field. 

The formalism of this Section is used in the next one 
to present the method of calculation of electroweak prop- 
erties of composite systems. Particularly, the method 
of construction of electroweak current operators is de- 
scribed. 



III. THE NEW RELATIVISTIC INSTANT-FORM 
APPROACH TO THE ELECTROWEAK 
STRUCTURE OF TWO BODY COMPOSITE 
SYSTEMS. 

In this Section we present our approach to elec- 
troweak properties of relativistic two-particle systems. 
To demonstrate how one describes the electromagnetic 
properties of composite systems in our version of the 
RHD instant form we first use the following simple model. 
We consider the system of two spinless particles in the S- 
state of relative motion, one particle having no charge. 
Let us note that a similar model was used in jlj where the 
authors gave the description of constituent interaction in 
IF of RHD and obtained the mass spectrum. The appli- 
cation of our method in general case follows the scheme 
of this Section. The case of tt - meson is investigated in 
Sec.IV and the S = 1 case in @. 

Electromagnetic properties of the system are deter- 
mined by the current operator matrix element. This ma- 
trix element is connected with the charge form factor 
F C (Q 2 ) as follows: 
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(pc \ump'c) = iPo+p%F c {Q 2 ) 



(26) reflections are: 



where p' c , p c are 4-momenta of the composite system in 
initial and final states, Q 2 = —t , q 2 = (p c —p' c ) 2 — t , q 2 
is the momentum-transfer square. The form (^6|) is de- 
fined by the Lorentz covariance and by the conservation 
law only and does not depend on the model for the in- 
ternal structure of the system. 

The Eq. (|2^) presents the simplest example of the 
extraction of a reduced matrix element. The 4- vector 
(Pc +Pc)ai describes symmetry and transformation prop- 
erties of the matrix element. The reduced matrix element 
(the form factor) contains all the dynamical information 
on the process described by the current. Usually, one 
does not fix the dependence of form factor on a scalar 
mass of the composite system p 2 — p' c 2 — M 2 , because 
it is diagonal with respect to this variable. The represen- 
tation of a matrix element in terms of form factors often 
is referred to as the parameterization of matrix element. 
The scattering cross section for elastic scattering of elec- 
trons by a composite system can be expressed in terms 
of charge form factor F C (Q 2 ). So, form factor can be ob- 
tained from experiment and it is interesting to calculate 
it in a theoretical approach. 

In this Section we calculate the form factor of our sim- 
ple composite system using the version of RHD IF based 
on the approach of the Section II. 

Now let us list the conditions for the operator of the 
conserved electromagnetic current to be fulfilled in rela- 
tivistic case (see, e.g., [fl0| ). 
(i). Lorentz- covariance: 



U- 1 (A)f(x)U(A)^Kr(A- 1 x). 



(27) 



Here A is the Lorentz-transformation matrix, U (A) - 
the operator of the unitary representation of the Lorentz 
group. 

(ii) . Invariance under translation: 



U- x {a)f{x)U{a)=f{x-a) 



(28) 



Here U (a) is the operator of the unitary representation 
of the translation group, 
(iii). Current conservation law: 



[&J"(0)]=0. 



(29) 



In terms of matrix elements ( j^(0) ) the conservation law 
can be written in the form: 



(30) 



Here q^ is 4- vector of the momentum transfer. 

(iv) . Current-operator transformations under space-time 



Up^fix^^Jix^x))^ 1 



U R j^(x)U^=^(-x) . 



(31) 



In (|31| ) Up is the unitary operator for the representation 
of space reflections and Up is the antiunitary operator of 
the representation of space-time reflections R = PT. 

(v) . Cluster separability condition: If the interaction is 
switched off then the current operator becomes equal to 
the sum of the operators of one-particle currents. 

(vi) .TVie charge is not renormalized by the interaction 
including: The electric charge of the system with in- 
teraction is equal to the sum of the constituent electric 
charges. 

In this paper the explicit equations for the form factors 
are obtained taking into account all the listed conditions. 



A. Electromagnetic properties of the system of free 
particles 

Let us consider first the simple two-particle system de- 
scribed in the beginning of Section III. The elastic scat- 
tering of a test particle, e.g., of an electron, by the system 
are defined by the operator of electromagnetic current 

jfi\o) of the two-particle free system. This operator 
can be calculated in the representation given by the ba- 
sis (||) or in the representation given by the basis @ . In 
the first case the operator has the form — ji^ ® 
Here is the electromagnetic current of the charged 
particle and I 2 is the unity operator in the Hilbcrt space 
of states of the uncharged particle. 



= (P2\p' 2 ){pi\3iM\p'i) 



(32) 



The matrix element of the one spinless particle current in 
the free case contains only one form factor - the charge 
form factor of the charged particle fi{Q 2 ): 



(pi\ji,(0)\p'i) = (pi+p'i),h(Q 2 ) 



(33) 



So, the electromagnetic properties of the system of two 
free particles ( p6| ) are defined by the form factor f\{Q 2 ), 
containing all the dynamical information on elastic pro- 
cesses described by the matrix element (3^) M. Particu- 
larly, the charge of the system is defined by the value of 
this form factor at Q 2 — > 0: 
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/i(0) 



(34) 



e c is the system charge. 

Now let us write the electromagnetic-current matrix 
element for the two-particle free system in the basis 
where the center-of-mass motion is separated (H): 



(35) 



(p,^,h; 0) (o) |P',V7) 



Here the variables which take zero values are omitted: 
J = S = I = 0. One can consider the matrix element j35| ) 
as a matrix element of an irreducible tensor operator on 
the Lorentz group and one can apply the Wigner-Eckart 
theorem. Under the condition of the theorem the matrix 
element of an irreducible tensor operator is the product of 
two factors: the invariant part (reduced matrix element) 
and the CG coefficient which defines the transformation 
properties of the matrix element. Thus, one can write 
(pq) in the form 



(36) 



(p,^\j [ ; } (o)\p\v^) = 

A ll {s,Q 2 ,s')g (s,Q 2 ,s') . 



The motivation for the parameterization ( pq ) is easy to 
be understood for our simple system. The 4-vector 
describes the transformation properties of the matrix ele- 
ment and the invariant function go(s, Q 2 , s') contains the 
dynamical information on the process. We will refer to 
go{s, Q 2 , s') as to free two-particle form factor. For more 
complicated systems the parameterization corresponding 
to the Wigner-Eckart theorem for the Lorentz group can 
be performed using a special mathematical techniques as 
described in the papers |52|, 34 1. 

The vector A ll (s, Q 2 , s') which describes the matrix- 
element transformation properties is defined by the 4- 
momenta of initial and final states only: we have no 
other vectors to our disposal. So A^(s,Q 2 , s') is a lin- 
ear combination of 4-momenta of initial and final states 
and is defined by the current transformation properties 
(the Lorentz-covariance and the conservation law): 



A„ 



s' + Q 2 )P^ + ( S '-s + Q 2 )P'}. (37) 



Thus, in the basis (Q) the electromagnetic properties of 
the free two-particle system are defined by the free two- 
particle form factor go(s, Q 2 , s'). 

So, in both representations (defined by the basis (||) as 
well as by the basis (||)) we pass from the description of 
the system in terms of matrix elements to that in terms 
of Lorentz-invariant form factors. 

One can see that (|||) and ( |36|) describe electromag- 
netic properties in terms of only one form factor. Both 



of these descriptions are, certainly, equivalent from the 
physical point of view. Let us consider the difference be- 
tween these descriptions. As we will show below by direct 
calculation the free two-particle form factor go(s, Q 2 , s') 
is not an ordinary function but has to be considered in 
the sense of distributions in variables s , s', generated 
by a locally integrable function. So, go(s, Q 2 , s') is a reg- 
ular generalized function. Let us remind that regular 
generalized function is that defined through an integral 
in the space of test functions. So, all the properties of 
go{s,Q 2 ,s') have to be considered as the properties of a 
functional given by the integral over the variables s , s' 
of the function go{s, Q 2 , s') multiplied by a test function. 
As test functions it is sufficient to take a large class of 
smooth functions that give the uniconvergence of the in- 
tegral. In particular, the limit ([m]) giving the total charge 
of the system through two-particle form factor is now the 
weak limit: 



Jim( 50 ( S , 



(38) 



Here </>(s, s') is a function from the space of test func- 
tions. The precise definition of the functional will be 
given below. 

As the invariant variables s , s' contain the energies of 
the relative motion of particles in initial and final states, 
one can consider the integral in ( j38| ) as an integral over 
these energies. 

At the first glance it seems that the description of 
the two-particle free system in terms of the form factor 
go(s,Q 2 ,s') is too complicated. However, so is the real- 
ity, as we will see later in the Subsection III.F. In fact, 
this kind of description is used implicitly for a long time 
in nonrelativistic theory of composite systems, without 
calling things by their proper names. It is this kind of 
description that makes it possible to construct the elec- 
tromagnetic current operator with correct transforma- 
tion properties for interacting systems. 



B. The form factor of the system of two free particle. 

The locally integrable function go(s, Q 2 , s') can be eas- 
ily obtained by use of CG decompositio n (jl9| ) for the 
Poincare group. Using (|l9|) we obtain for (|36|): 



dpx dp 2 dpi ' dp 2 ' , - i — 1 _ . 



2p w 2p 2 o 2p' 10 2p' 20 
x (pilP2\j^HQ)\p[;p' 2 )(p[-,f 2 \P', v^ 7 } • (39) 
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To calculate the free two-particle form factor one has to 
use ([32k (|33|), ( |36| ) and the explicit form of CG coeffi- 
cients ( 19j) for quantum numbers of the system. As the 
particles of the system under consideration are spinless, 
now ([l9]) does not contain D - functions. 

It is convenient to integrate in ( |39| ) using the coordi- 
nate frame with P' = , P = (0, 0,P). As the result 
we obtain the following relativistic invariant form for the 
function go(s, Q 2 ■, s'): 

9o(s,Q 2 ,s') 
(s + s' + Q 2 ) 2 Q 2 $(s,Q 2 ,s') 



2V(s-4M2)( s '-4M2) [A(s,-Q2,s')] 3/2 



h(Q 2 



(40) 



Here tf{s,Q 2 ,s') = 6{s' - s x ) - 0(s' - s 2 ), and is the 
step function. The result, naturally, does not depend on 
the choice of the coordinate frame. 



si, 2 = 2M 2 + ^p(2M 2 + Q 2 )(s - 2M 2 ) 



2Af 



:^Q 2 {Q 2 + 4M 2 )s{s - AM 2 ) 



The functions s\.2{s,Q 2 ) give in the plane (s, s') the 
kinematically available region. The position of this re- 
gion depends strongly on the momentum-transfer square 
t = —Q 2 0. The simplest way to obtain the functions Si^ 
is a geometrical one [B3|. 




(p-n 

FIG. 1. Kinematical triangle 



2 The $ - function is a purely kinematical factor for IA. 
This fact does not depend on relativism and takes place in 
nonrelativistic case, too. See Subsection III.F for details. 



Let us construct the triangle schematically presented 
on the Fig.l. The side OB presents the vector pa^, and 
p' 2 , as following ( |32] ) p 2fl =i>2u- The side CB presents 
the vector pi^ and AB - p' x . Now the sides of the 
large triangle AOC present the vectors P^ , P M and (P M — 

P'), with the norms Vs 7 , \J~s-, Vi, respectively. Vectors 
of initial state p[ , p%^ = p' 2fl , P^ are fixed. Let 
us fix the norm of vector (P M — P' ). So , because the 
norm of vector pi^ , p\ = M 2 is constant, the triangles 
ABO and ABC are determined unambiguously by three 
sides. But the triangle ABC can be rotated around the 
side AB (p' lAt ). It is possible to find the minimal ^fsi 



and maximal y/s2 lengths of OC (norm of vector P^) 
under this rotation. The value of the Si , s 2 give the 
kinematically available region in the plane (s , s') which 
is symmetric under interchange s <-» s'. 




FIG. 2. The kinematically available region in the plane s, s' 
(inside the parabolae). The calculation is performed for: 1) 
Q 2 = 2M 2 . 2) Q 2 = Af 2 /2. 3) Q 2 = Af 2 /64. The constituent 
mass is M = 0.25 GeV. 



In Fig. 2 the domain where the generalized function 
( plO| ) is nonzero in the plane s , s' is given for different 
values of the momentum-transfer square. One can see 
that the free two-particle form factor g (s,Q 2 ,s') ( |4fj| ) 
has in fact to be interpreted in terms of the distributions: 
The ordinary limit as Q 2 — > is zero because of the 
cutting i? - functions and the static limit exists only as 
the weak limit (j38|) . 

Let us calculate this limit. Let us define the functional 
giving regular generalized function as a functional in R 2 
as follows: 
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(<7oO, Q 2 , s') , (p(s, s')> = J dfj,(s, s') g (s, Q 2 , s') <p(s, s') . 

(41) 

Here 

g (s,Q 2 ,s') 

= 16 6{s - AM 2 ) 6{s' - AM 2 ) g (s, Q 2 , s') . (42) 



<i/i(s, s') = wsP cfyx(s) dyu(s') 



dfj,(s) = ifcdVs • (43) 



The 9 - functions in (|4^) give the physical region of possi- 
ble variations of the invariant mass squares in the initial 
and final states explicitly. The measure (^) is due to 
the relativistic density of states (p2|), (p5|). cf>(s,s') is a 
function from the test function space. So, for example, 
the limit of go(s , Q 2 , s') as Q 2 — ► (the static limit) has 
the meaning only as the weak limit (compare with (p4j)): 



]im o (ff ,^> 



= (eS( M (s') - n(s))$(s - AM 2 ),4>) . (44) 



It is this weak limit that gives the electric charge of the 
free two-particle system. If the test functions are nor- 
malized with the relativistic density of states, then the 
r.h.s. of the Eq. (Q) is equal to the total charge of the 
system. 



C. Electromagnetic structure of the system of two 
interacting particles. 

Now let us consider the electromagnetic structure of 
our simple model (^6|) in the case of interacting particles. 

As we have mentioned in Sec. II when constructing the 
bases (JsJ) and (^) in the frame of RHD the state vector 
| p c ) belongs to the direct product of two one-particle 
spaces. We can write the decomposition of this vector 
with J = 1 = S = mj =0 in the basis (|). Now @ has 
the form: 



dPrIP' _ _ 

" —dV~sd^( Pc \p,V^)(P,V^\Uo)\P',^) 



NcgN>c G 

x (?' , sfi>\p' c ) = {Pc+p' c )»F c {Q 2 ) . (45) 



Here (P ' , v s'\p' c ) is the wave function in the sense of the 
instant form of RHD @. 
Using (|l|) we obtain for (||): 

J McG J^CG 



x (pc , V^U^OM ' , Vs 7 ) = (p c + Fc(Q 2 ) ■ (46) 

We have omitted in the wave function^ 23) the variables 
with zero values: J = S = I =0 (se c (p5 ) too). 

Using ( 22|) , (O) we can rewrite (f46|) in the form of the 
functiona in Rr: 

dn{s , a') u(k(s)) J M (p c , \/s\p c ' , Vs') u(k(s')) 

= (Pc+p' c ) t ,F c (Q 2 ) , (47) 
J„(p c ,Vi;p c ',^) = 16 6(s-AM 2 )9(s' -AM 2 ) 



N C N> 
N CG N' CG 



(Pc,Vs\jn\Pc' , 



In the previous cases the state vectors and the op- 
erators entering matrix elements transformed following 
one and the same representation of the Poincare group 
(namely, following the universal covering subgroup of the 
Poincare group - nonuniform group SL(2, C) Now 
in the matrix element in the integrand of ( p!7| ) the state 
vectors and the operator transform following the differ- 
ent representations of this group. The current operator 
describes the transitions in the system of two interact- 
ing particles and transforms following the representation 
with the generators of Lorentz boosts depending on the 
interaction (0). The state vectors belong to the basis 
(^) and physically describe the system of two free par- 
ticles and, so, transform following a representation with 
generators which do not depend on the interaction (|2C|) . 
That is why the current operator matrix element in ( |47j ) 
can not be represented in the form (p^), (pTl): we can 
not construct the 4-vector defining the matrix-element 
transformation properties under Lorentz boosts from the 
variables which the state vectors depend on. 

Nevertheless, as we show below, the problem of the 
parameterization of the current matrix element in ( |47| ) 
can be solved if one consider this equality as the equality 
of two functionals. The l.h.s. contains a functional in 
R 2 generated by the Lorentz-covariant function (current 
matrix element). Let us denote 



i/>(8, s')=u(k{s))u(k'(s')) 



(48) 



The functional in the l.h.s. of ( p7| ) is given on the set of 
test functions ip(s , s') through an integral in R 2 and de- 
fines a Lorentz-covariant (regular) generalized function 
with the values in the Minkowski space (see, e.g., El). 
Here Q 2 is a parameter. The test-function space can be 
(in general) larger than (^). However, the uniconver- 
gence of ( ]47| ) has to be guaranteed. 
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Let us write the matrix element in the form analogous 



Wc^Pc'.V?) = B„{s,Q 2 ,s')G(s,Q 2 ,s') . (49) 

The covariant part in (]4^) (as well as in (|36|)), the vec- 
tor Bfj,(s , Q 2 , s'), is supposed to be an ordinary smooth 
function and the invariant part G(s , Q 2 , s') is general- 
ized function. In fact, G(s , Q 2 , s') is the reduced ma- 
trix element containing the information on the process. 
This kind of representation of a Lorentz covariant gener- 
alized function as a product of a Lorentz covariant ordi- 
nary smooth function and a Lorentz invariant generalized 
function was described in jyj. Let us remark that all 
normalization constants enter the vector B^. Contrary 
to (|36|), now it is impossible to construct the 4-vector 
B^, describing the transformation properties of the ma- 
trix element, in terms of the variables entering the state 
vectors (as was pointed out before). 

Using (E9J) we can rewrite (E7J) in the following form: 



d»(s , s') V(s , s')B M (s , Q 2 , s')G(s , Q 2 , s') 

= (pc+p' c )^F c mQ 2 ) . 



(50) 



To obtain the vector B^ let us require the Eq.(^0]) to be 
covariant in the sense of distributions, that is to be valid 
for any test function ip(s , s') in any fixed frame. The 
variation of test function in the functional (50) means in 
fact, following (|^|), the variation of the wave function of 
the internal motion. Under such a variation the vector 
in the r.h.s of ( |50| ) is unchanged as it is constructed with 
4-vectors describing the motion of the system as a whole, 
independent of the internal constituent motion. As to the 
form factor in the r.h.s it varies under the test function 
variation. So, under a variation of the test function the 
r.h.s. of ( pp[ ) remains to be collinear to the vector (p c + 
p' c )fj,. At the same time, under arbitrary variation of the 
test function the vector in the l.h.s. in general changes 
the direction. So, for the validity of the equality ( |50| ) 
with arbitrary test function it is sufficient to require that 
the following equation 



B^s,Q 2 , S ') = ( Pc +p' c ), 



(51) 



holds. This choice of the vector B^ in (|5l|) ensures that 
the l.h.s. of ( |47| ) satisfies the condition of Lorentz covari- 
ance for the current as well as the condition of current 
conservation. 

Let us discuss the physical meaning of the representa- 
tion (|49|), ( fill ) for the matrix element. As this represen- 
tation is explicitly Lorentz covariant and also satisfies the 
current conservation law, then it means that the current 



operator for the composite system contains the contribu- 
tion not only of one-particle currents but of two-particle 
currents, too (see, e.g., 0)): 



i(z) = £i (fe) (z) + £ 



(ki 



(52) 



Here the first term is the sum of one-particle currents and 
the second - of two-particle currents. In the case of our 
simple model each sum in ([52] ) contains only one term. It 
is well known that if one approximates j(x) « ^2 k (x) 
then the current operator in IF dynamics does not satisfy 
the condition of Lorentz covariance and the conservation 
law So, from the physical point of view, the covariant 
part of the current matrix clement (51) which defines the 
transformation properties of the current in ( f47| ) is given 
by (^2|) and contains the contributions of one- and two- 
particle currents. 

The invariant part of the decomposition (f49|) is the 



form factor or the reduced matrix element Gl 



and contains the information on the dynamics of the scat- 
tering of test particle by each of the constituents (the first 
term in (|52"|)), i.e. by the free two-particle system, as well 
as by two constituent simultaneously (the second term). 
So, the form factor contains the contribution of the free 
system form factor ( p0| ) and the contribution of some ex- 
change currents analogous to meson currents in nucleon 
systems M. 



G(s , Q 2 , s') = g (s, Q 2 , s') + G c (s , Q 2 , s') 



(53) 



Here G c is the reduced matrix element containing the 
contribution of two-particle currents ([52]). 

Using (||), (||), (|8l), (|l|) one can obtain from @ 
the scalar equation of the following form: 



dV^dV?v(s)G(s,Q 2 ,s')v(s') = F C (Q 2 ) . (54) 



The form factor G(s , Q 2 , s') includes all possible mech- 
anisms of the transition described by the matrix clement 
(ffif). So, the representation (^i|) for the charge form fac- 
tor of the system is quite general. 

Now let us proceed with the approximate calculation 
of the form factor (|54|). 



D. Modified impulse approximation (MIA) 



The problem of the calculation of the form factor 
G(s , Q 2 , s') (|54|) including exchange currents is a very 
difficult problem. We propose an approximation which 
is a kind of analog of relativistic impulse approximation. 
We propose to omit the contribution of the two-particle 
currents to the form factor G(s , Q 2 , s 1 ). 
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However we will not change the covariant part S M of 
the current matrix element in (fi"9|), so that this covari- 
ant part will contain the contribution of the two-particle 
currents and so that the transformation properties of the 
matrix element will not be changed. 

So, we approximatel y c hange the generalized function 
G(s , Q 2 , s') in ( [49|), ( p3|) for the generalized function 
go(s, Q 2 , s') (p6|), (|40[) , which describes, as we have shown 
before, the electromagnetic properties of the free two- 
particle system. Nevertheless, the matrix element (45), 
( [49| ) as a whole will contain the contributions of two- 
particle currents, although not the full contribution but 
such that ensures its correct transformation properties. 

Let us note that our approximation does not contra- 
dict general statements (see 0) that to obtain correct 
description of electromagnetic current of composite sys- 
tem which satisfy the Lorentz-covariance condition and 
the current conservation law one has to take into account 
many-particle currents. 

Let us discuss now the meaning of our approximation 
from the point of view of the Wigner-Eckart theorem 
for the Lorentz group. The matrix element of a current 
including many-particle currents, following the Wigner- 
Eckart theorem for the Lorentz group, can be presented 
in the form (ff9|), (|5l|). The dynamical information 
on many-particle currents is contained in the reduced 
matrix element - the form factor, while the transforma- 
tion properties of the contributions of many-particle cur- 
rents are defined by the covariant part of the form (49). 
So, our approximation means that the dynamical part 
of the contribution of the many-particle currents to the 
total current is omitted while the covariant part of the 
contributions remains. The dynamics of many-particle 
currents remains out of the limits of the approximation, 
while the transformation properties of the total current 
remain intact. 

Thus, in our approximation the scalar equality ( p34| ) 
transforms into approximate scalar equality which cor- 
responds, from the physical point of view, to relativis- 
tic impulse approximation. In the developed mathemat- 
ical formalism we have not broke the Lorentz covariance 
of the current nor the current conservation law. Let us 
point out that to calculate form factor we do not use a 
special current component as it is done in other mathe- 
matical formulations of RHD (see, e.g., Q). Let us re- 
mark that, from the physical point of view, the form fac- 
tor go(s, Q 2 , s') contains the contributions of one-particle 



F(Q 2 )= / d^dV?<p(s)go(s,Q 2 ,s')v( S ') . (55) 



currents only (see Equations (|36|), (|39|), fllpP) and in this 
sense our approximation corresponds to the known im- 
pulse approximation. In order to emphasize that our ap- 
proximation differs from the usual IA we will refer to it 
as to modified impulse approximation (MIA). The form 
factor of the composite system in MIA has the form: 



It is worth to notice that the Eq. (|5l|) and the form ( pq ) 
can be formally obtained if we write in ( [45| ) the current 
of the free system ( |36| ) instead of that of the interaction 
system and change the covariant part of (|37]) for 



(56) 



A f _ l (s,Q 2 ,s / )\ P=Pc , P i= p , c = {p c +p' c )^ 



The Eq. (|56|) gives a simple prescription to write the 
current matrix elements for interacting system in the ba- 
sis (j^) in MIA using the current parameterization ( |36| ) 
for the free system. The prescription is as follows: in 
the vectors in the parameterization (|36|), ( j37|) one has 
to use the momenta of composite system instead of the 
center-of-mass momenta of the free two-particle system. 
Note that this prescription works for more complicated 
systems, too. 

We do not discuss in this paper the problem of going 
beyond the limits of MIA and of obtaining corrections to 
go {s , Q 2 , s') in (l53l), (|55|). This means that if considering, 
for example, nucleon systems we do not take into account 
meson current. 

Let us consider now the fulfilling of the conditions (i)- 
(vi) for the electromagnetic current. 

The conditions (i)-(iii) are satisfied by construction. 
For example the fulfilling of (i) and (iii) is ensured by 
the correct transformation properties of the 4-vectors in 
@, ©, and ©. 

The condition (iv) is satisfied immediately as the 
form factor <7o(s , Q 2 , s') in ( |36"| ) and the form factor 
G(s , Q 2 , s') in ([4"^ ) are scalars in our simple model 

The condition of cluster separability (v) needs a more 
detailed consideration. At large distances (or if the 
interaction is switched off) the contribution of two- 
particle currents has to go to zero: G c (s ,Q 2 ,s') — > 
in (p3|). This means that in the form (^3h the form fac- 
tor G(s , Q 2 , s') has to transform into go(s , Q 2 , s'). Let 
us remark that the condition of cluster separability is 
fulfilled in MIA, too, as in this approximation the use 
of go(s , Q 2 , s') instead of G(s , Q 2 , s') is supposed from 
the very beginning. When the interaction is switched off 
the generalized function g {s,Q 2 ,s') for the free two- 
particle system acts on a larger space of test functions 
than (|48|). As go(s , Q 2 , s') contains only the one-particle 
current contributions (p9) the condition (v) is satisfied 



3 The currents which do not conserve the parity also can be 
considered in our formalism. In that case one can separate 
not only the scalar part of current matrix element but the 
pseudoscalar part, too. This case is considered elsewhere. 
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and the composite-system current go to the sum of the 
one-particle currents. 

The condition on the charge to be nonrenormalizable 
also is fulfilled directly in MIA because the weak limit 
( [lif ) does exist on test functions (ffil). 

So, our prescription for the construction of the cur- 
rent in MIA satisfies all the conditions for the current 
operator. 



E. MIA versus IA 

Let us compare the approximation MIA with the well 
known IA. 

To do this let us first calculate the form factor in IF 
RHD not using the canonical parameterization. In par- 
ticular, let us formulate IA in terms of operators as it is 
formulated usually (not in terms of form factors). Let us 
decompose the matrix element ( p6| ) through the complete 
set of states (0): 



P. 



Um(o)Ipc> = 



dpi dp 2 dpi ' dp 2 ' 
2pio 2p 2 o Vio 2^2 



2H 



(pVspVb'c) 



(57) 



Here ( p\; p% \p c ) is wave function of constituents in com- 
posite system. 

If the current matrix element in (j57j) is taken in the 
IA approximation ( |52"| ) and contains one-particle cur- 
rents only, then the Eq. ( ^7|) is selfcontradicting Q . In 
fact, one can show this in the following way. In our sim- 
ple model all the dynamical information about the cur- 
rent (i.e. the composite system form factor) can be ob- 
tained from only one matrix element in the Brcit frame. 
However, to go to the Breit frame one has to perform 
the transformation which is interaction depending. This 
means that in Breit system two-particle currents appear 
along with one-particle ones. The form factors calculated 
in arbitrary coordinate frames using different matrix el- 
ements will be of different forms. 

To write the form factor in terms of wave functions 
( pl| ) one has to perform the CG decomposition of the 
basis (||) in terms of the basis (^|) in the wave functions 
( p7| ) and to use the explicit form for CG coefficients (|H 
for the quantum numbers of the system: 



[pr,P2 \p c ) = \l^{P,Vs,J,l,S,mj \p c 



(58) 



The current matrix element in ( |57| ) has the form (32) 
The one-particle currents are expressed through the form 
factors (I3 - 



The Eq.(p7[) is an equality for two 4- vectors. Taking 
different components of this equality and exploiting 5— 
functions in integrals, one can calculate the form factor 
of the composite system. The result of calculation of the 
form factor in this way is not unambiguous. In particular, 
it depends on the actual choice of the component of the 
current (|5^) to be used in the calculation. Moreover, the 
result depends on the coordinate frame chosen to perform 
the integration in (57). This is the general feature of IA 
in the usual formulation of IF RHD (see, e.g., 0). 

Let us write the final result of the calculation of the 
form factor from the equation for the null-component of 
the current and performing the integration in the coordi- 
nate frame where p c ' = , p c = (0, 0, p). If now we write 
the integral in terms of the invariant variables s, s' the 
obtained form factor has the form: 



F C (Q 2 ) 



M c ^2(2M 2 + Q 2 ) 



4M 2 + Q 2 



d\fs d 



(s + s' + Q 2 ) 4 



s' v /(s-4M 2 )(s'-4M 2 ) [A(s , - Q 2 , s')] 3/2 

i e{ Sl Q 2 ,s'^ 



( SS ') 1/4 \/^(s + Q 2 ) 



^(s)^( S ')/i(Q 2 ) ■ (59) 



The Eq.(|59|) differs from (p5|), obtained with the use of 
the two-particle free form factor. In the case of wave 
functions satisfying the conditions j2^), (|23|), the form 
factor ( |59| ) satisfies the normalization: F c (0) = e c . Let 
us note that the form factor obtained in this way from 
the third current component in ( [j7| ) does not satisfy this 
condition. 

Let us compare IA and MIA results. Let us note once 
again that in MIA we separate (by use of the scheme of 
canonical parameterization) the covariant part of the cur- 
rent matrix element in (pO) prior to perform any calcu- 
lations. This covariant part ensures the correct transfor- 
mation properties of the corresponding decompositions 
in terms of free-particle states. The difference between 
@ and @ is: 



xg ( S ,Q 2 , S ') [l-R(s,Q 2 , S ')] 



(60) 



R{s,Q 2 ,s') = 



M c v/2(2A/ 2 + Q 2 ) 



2 4M C 2 + Q 2 

(s + s' + Q 2 ) 2 1 

(ss') 1/4 V^(s + Q 2 ) 



(61) 
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The value R(s, Q 2 , s') presents an additional factor to 
one-particle currents, that is in reality the two-particle 
current contributions. This term ensures the Lorentz co- 
variance of the electromagnetic current matrix element 
and the current conservation law in (|4^ ) . Let us note that 
this additional term contains no dynamical information 
on the interaction of test particle with two constituents 
simultaneously. It does not depend, for example, on the 
interaction constants for such a process. 

So, to summarize, we can write the following schematic 
equations: 

(IA) Brett 



(MIA) BreU = (MIA) Lab 

It is well known that the standard IA depends strongly on 
the coordinate frame used for the calculation. The MIA 
results do not depend on it at all. So, the differences 
between IA and MIA results for different IA coordinate 
frame can be rather significant. 

Notice that IA and MIA coincide in the nonrelativis- 
tic limit. As this takes place, the nonrelativistic limits 
of form factors, which were obtained from the different 
current components, are identical. Hence the difference 
between IA and MIA is connected with the breaking of 
relativistic covariance conditions really. 

We give the quantitative comparison of the form fac- 
tors obtained in IA and in MIA in the Section IV where 
the realistic calculation of the pion electromagnetic struc- 
ture is given. 



F. The nonrelativistic limit 

The description of composite-system form factors in 
terms of distributions is not a specific feature of our rel- 
ativistic approach. The similar formalism is used in non- 
relativistic theory of composite systems |12| for a rather 
long time (although not referring to the mathematics of 
distributions). In the nonrelativistic limit our approach 
gives the formalism developed in j42] . 

In the nonrelativistic limit the relativistic charge form 
factor (153) has the following form: 



F N r(Q 2 )= / k 2 dkk /2 dk'u(k)g 0NR (k,Q 2 ,k')u(k') 



gmR{k,Q\k') = f -0^e{k, 



X) 



(62) 



(63) 



${k,Q 2 ,k') = 9[k' 



Q 
2 



9 k' 



Here goNn{k,Q 2 , k') is the free relativistic form factor 
obtained from ( pp| ) in the nonrelativistic limit. fi(Q 2 ) 
is the charged-particle form factor. The obtained result 
coincides with that derived in standard nonrelativistic 
calculations E). 

In p2] | the same formulae are obtained from the equa- 
tions for form factors in terms of coordinate representa- 
tion wave functions: 



f N r(Q 2 ) = h(Q 2 



drr ir{r)jn ( ^ 



(64) 



The Eqs.(p2[), ([33]) can be obtained from (64) by use of 
the Bcsscl transformation: 



u(k) 



r dr u(r) jo (kr) 



(65) 



and the normalization condition: 



u (r) dr 



k 2 u z {k) dk = 1 



Rigorously speaking, the Eq. (|62j ) has to be interpreted 
as a functional in the sense of distributions generated by 
the function goNH.(k, Q 2 , k') and defined on test functions 
u(k) u(k'). The ordinary function ( p3[ ) generates regular 
generalized function defined generally on the larger class 
of test functions ip(k,k') in R 2 , providing the uniform 
convergence of the integral. One needs the uniform con- 
vergence to take limits in the integrands. 

Let us define the functional in R 2 by the following 
regular distribution (compare with (41)-(p3)): 



{goNR(k,Q 2 ,k'),^(k,k')) 
dfi(k, k') g 0NR (k, Q 2 , k') 1>(k, k') , (66) 



g N R (k, Q 2 , k') = tf(fc) #(k')g mR (k, Q 2 , k') , 

d/i(k, k') = d/i(k) dfj,(k') , d/i(k) = k 2 dk . 

The function goNR{k,Q 2 ,k') which appears in H|] 
quite formally, here has a definite physical meaning 
and describes the electromagnetic properties of non- 
relativistic free system of two spinless particles in the 
S - state, one of particle having no charge (compare 
with g (s,Q 2 ,s') in d|), @, @). The static limit 
lmiQ2^ goNR{k,Q 2 ,k') giving the system charge exists 
only in the weak sense as the limit of the functional (|6^) : 

lim {goNR(k,Q 2 ,k') ,<ip(k,k')) 
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(67) 



On the test functions ip(k, k') = u(k) u(k') (with u(k) 
- being the normalized bound state wave function), the 
functional ( |66| ) defines the bound state form factor in the 
nonrelativistic IA (|6^). The weak limit (|67]) is equal to 
the system charge: 



lim (g NR(k,< 
Q 2 — >0 



\k')^(k,k')) 



k 2 dk u 2 {k) 



(68) 



So, one can see that the description of the system in 
terms of form factors in IA by the Eq.(|62|) (as in filfl ) in 
fact defines the form factor in the sense of distributions 
as a functional defined on the set of the wave functions 
of bound system. 

To go beyond nonrelativistic IA one has to addend 
some terms to goNFt(k, Q 2 , k'). For example, such terms 
cause the meson exchange currents in two-nucleon sys- 
tems. So, in the standard nonrelativistic theory the dy- 
namical treatment of exchange currents is performed in 
the same way as in our relativistic approach (p3|). 

So, to conclude, one can consider our approach to IA to 
be a relativistic generalization of nonrelativistic IA, and 
our equations for form factors in this approximation to be 
a relativistic generalization of the equations of {42| . Let 
us remark that in more complicated systems (e.g., in p 
- meson and deuteron) our relativistic form factors also 
have correct nonrelativistic limits which coincide with 
0. 



G. A bridge to dispersion relations 

Let us discuss now one of the unsolved problems of 
RHD - the possible links between RHD and quantum 
field theory (QFT) f§. The fact, that RHD, contrary to 
QFT itself, operates with the finite number of degrees of 
freedom, makes it to be in some way similar to the disper- 
sion approach of QFT, which is dealing in principle with 
a finite number of degrees of freedom, too. However, the 
dispersion relations, based on the analytic properties of 
the scattering amplitudes, matrix elements, form factors 
in the complex energ y p lane, are rather correctly derived 
in the frame of QFT [43]. So, it seems to us, that one can 
look for links between RHD and QFT not only directly 
but through the dispersion approach, too. 

Here, using the simple model of the previous Subsec- 
tions, we compare our version of RHD with the so called 
modified dispersion approach. Dispersion-relation inte- 
grals over composite-particle mass are used in this ap- 
proach. This approach enabled one to write the deuteron 



electromagnetic form factors in terms of the physical 
hadron scattering phase shift and gave the results for 
the elastic ed-scattering in good agreement with experi- 
mental data. The details of the modified dispersion ap- 
proach can be found in p3| , ^ -^6| (see also jl7)). Let us 
note that an immediate application of the approach to 
quark systems is difficult to realize because of the fact 
of quark confinement. However, there are some inves- 
tigations based on similar ideas where the form factors 
of hadrons as constituent-quark bound states are con- 
sidered in the frame of the dispersion technique of the 
integral over composite particle mass p8j. 

For convenience of reader let us describe briefly, omit- 
ting the proofs, the essence of the modified dispersion 
approach and obtain the electromagnetic form factor of 
composite system for our simple model following the pa- 
per ||. 

The Heisenberg current operator of the system of two 
particles interacting as in (p3) can be written in the form 



- ,-(0) , Aint) 



(69) 



Here is the current of the free two-particle system 



(see (32), (|36|)). The operator j^ nt ^ is interaction depen- 
dent. Let us suppose that our model constituent system 
has scattering states and let us calculate the matrix ele- 
ment of the operator ( |69| ) between in- and out- states. 
Let us suppose that the scattering states are the S- states 
of relative motion. The matrix element of the operator 
jjx can be written in terms of the free two-particle form 
factor (|3^) calculated previously (flO|). The matrix ele- 
ment of the operator containing the interaction 



(P(±) I ii int) 



p'(±)) 



can be written in terms of the form factor 
Gi(s =F is, Q 2 , s' ± ie) . 



(70) 



Here sign "+" stands for in- state and sign "— " for out- 
state. The form factor in (|7C| ) has kinematic cuts in the 
complex plane of the variables s , s' . The cuts go along 
the real axis from the point 4 M 2 to infinity. The notation 
Gi(s + ie, Q 2 , s' — ie) means that there exists the analytic 
continuation of this form factor from the physical region 
of the variable s into the upper complex half-plane, and 
to the lower half-plane in s' . One can check this fact 
considering simple models. The form factor entering the 
parameterization of the total current ( |69|) can be written 
as the sum of the form factors: 

G(s, Q 2 ,s') = go( S) Q 2 ,s') + G,( S , Q 2 , a') . (71) 

Let us consider the matrix element of the total current. 
Let us fix the variable s and let us connect by i'-matrix 
the in- and out- vectors of the basis in the variable s': 
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(P\3»\P'(+)) = {P\j»\P'(-))S( S i). (72) 

S(s) = cxp(2i<5), S is the scattering phase shift. S - 
matrix can be written in the form: 



S(s) = 



B(s-ie) 
B(s + ie) 



(73) 



Here B(s) is the relativistic analog of the Jost function. 
Taking into account Eqs . (ffoj) , ( [73] ) one can rewrite fl7^ ) 
as: 

Gi{s, Q 2 ,s' - ie) B{s' - ie) - G;(s, Q 2 ,s' + ie) B(s' + ie) 
= - g (s, Q\ s') (B(s' - ie) - B(s' + ie)) . (74) 



The equation ([74]) presents the so called Riemann- 
Hilbert problem for half-axis. The solution has the form 
j§: 

Gi{s, Q 2 , s') B(s') = G{s, Q 2 , s') + d{s, Q 2 ,s') , (75) 



G(s,Q 2 ,s') = - 



ds"g (s,Q 2 ,s")A(s") 



2ni 



'AM 2 



(76) 



A(s) = {B(s + ie) - B(s - ie)) , 

Here Ci(s,Q 2 ,s') is an unknown function, regular in s' 
in the neighbourhood of the real axis for AM 2 < s' < oo. 
Now let us connect the in - and out - bases in the variable 
s. Taking into account the explicit form of Gi (75), ( |76| ) 
we obtain the boundary value Riemann-Hilbert problem 
for the function C\ in the variable s with the solution: 

Ci (s - ie, Q 2 , s') B(s - ie) - C x (s + ie, Q 2 , s') B(s + ie) 



= g (s,Q 2 ,s')B( s ')+G( S ,Q 2 , S ') A(s) 



(77) 



B(s) d( S ,Q 2 , S ') 



2ni 



ds" 



AM 2 



g {s"',Q 2 ,s')B{ s ') + G{s'" 1 Q 2 lS ') A(s'») 



s — s 

2 „/ 



/// 



+ C{s,Q\s') 



(78) 



The unknown function C(s, Q 2 , s') is regular in the neigh- 
bourhood of the real axis for AM 2 < s,s' < oo in s' as 
well as in s. Now let us consider the matrix element of 



the total current 
ing @, ©, ([ 



(1691) , (yW between in-states. Substitut- 
in ( [71[ ) we obtain finally the following 
form for the form factor of the total current (|6^) in in- 
basis: 



G(s,Q 2 , S ')=g (s,Q 2 ,s') 



1 



2niB(s' + ie) J 4 



ds 



,>90{S, 



,s")A(s") 



M 2 



ie 



1 



2 ni B(s — ie) 



ds' 



,g (s>",Q 2 , S ')A( s >") 



AM 2 



X / ds' 1 
Um 2 



Air 2 B(s-ie)B(s' + ie) 



°° dsll go( S "',Q 2 , S ")A( S "')A(s") 



4m= (s-s'»-is)(s> 

C(s,Q 2 ,s') 
B(s-ie)B{s' + ie) ' 



s" + i e) 



(79) 



The Eq.(f79|) provides the correct analytic properties of 
the form factor obtained in QFT approach J43J . In partic- 
ular, this form contains the anomal branch points known 
from the dispersion approach to composite systems (e.g., 
to the deuteron). The Eq.([79|) can be used to obtain the 
form factor of the constituent bound state for the case 
of the 5-state of relative motion. Now it is necessary 
to perform the analytic continuation of ( f79| ) in the vari- 
ables s,s' to the bound state point s = s' = M 2 (M c 
is the bound state mass) and to take the residues in the 
poles. As the result we obtain the bound-system form 
factor directly in terms of the S'-scattering phase shift 
for constituents: 



F C (Q 2 



r- 



ds ds 



AM 2 



,g (s,Q 2 ,s')A(s)A(s') 
(s-M 2 )(s>-M 2 ) 



(80) 



Now the constant T 2 is determined by the condition 
-Fc(O) = e c and indirectly takes into account the contribu- 
tions of the so called unphysical cuts. The Jost-function 
discontinuities can be written in terms of experimental 
scattering phase shift. The free two-particle form factor 
for our model has the form (40). F C (Q 2 ) is the functional 
generated by the generalized function <?o(s, Q 2 , s') on the 
test functions A(s)/(s — M 2 ). The described formalism 
was applied to the deuteron in Ref. jl6| and gave a good 
agreement with experimental data. 

Let us note, that the form ( |80| ) obtained through the 
modified dispersion approach is in close analogy to the 
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forms @, @ obtained in the frame of IF RHD. This 
analogy can be made even more obvious using the results 
of the Ref. [^5| where neutron-proton system was con- 
sidered in nonrclativistic case. In (4^] it is shown that 
if the deuteron electrodisintegration amplitude satisfies 
Mandelstam representation, then the wave function of 
the system has the well fixed form and can be expressed 
in terms of np -scattering phases. This wave function 
satisfies a dispersion relation. The analytical properties 
used during the derivation of this relation are the same 
for the large group of phenomenological potentials. So 
the obtained dispersion relation can be used to find ex- 
plicit form for real two-nucleon systems wave functions. 
As one can see from the solution structure, such a re- 
construction of wave function is stable both in the usual 
sense and in the sense of the large energy phases influ- 
ence. Finally, the bound state wave function u(r) is of 
the form: 

u(r)=f / dx A( - X \ [sm(xr)}. (81) 

J-oo [X-K) 

(See for details.) Here k 2 is the deuteron binding 
energy and the nonphysical cuts contribution enters the 
normalization constant T. 

In nonrclativistic case the Jost-function discontinuity 
multiplied by the pole term in (^o|) gives the nonrelativis- 
tic wave function up to nonphysical cuts contribution 
fj5| . When these contributions are taken into account 
the equations (|2|) and ( [SOl) do coincide. 

So, the relativistic Eq.(|55|) can be motivated (at least 
for composite systems which have the scattering states) 
in the frame of modified dispersion approach, that is on 
the usual level of correctness for obtaining the analytical 
properties of the form factors of composite systems in the 
frame of QFT. 

IV. THE ELECTROWEAK STRUCTURE OF 
PION 

Now we apply the method of previous sections to the 
calculation of the electroweak structure of pion. There 
exists a lot of experimental data on pion, so the effec- 
tiveness of the method can be checked by the comparison 
with the data (see, e.g., Jl6| and references therein). 

A. The electromagnetic form factor of pion 

There are many facts that make it interesting to con- 
sider the pion in the frame of the formalism developed 
in the previous sections. First, the pion is an important 



object in the particle physics and is in the focus of inter- 
est for years. Second, the pion consists of light quarks 
and thus has to be considered in the frame of relativistic 
approach. Third, at the present time the large program 
on the pion structure is on line in JLab ]5(J. 

We consider the pion as a system of two constituent 
quarks. The system is described by a phenomenological 
wave function. 

In theoretical treatment of the pion electromagnetic 
structure one has to make difference between "soft" and 
"hard" parts of the form factor. The "soft" part which 
dominates at small and intermediate momentum trans- 
fers, needs nonperturbative approaches. The " hard" part 
which defines the form factor at asymptotically large val- 
ues of momentum transfers can be calculated from per- 
turbative QCD. However, a controversy still exists con- 
cerning the scale of momentum transfers characteristic of 
the transition from the nonperturbative to the perturba- 
tive regime (see, e.g. [[[l]). It is pointed out by different 
authors that the existing experimental data are defined 
by the "soft" part of the form factor and can be described 
by use of phenomenological wave functions of constituent 
quarks without involving the perturbative QCD. Usually 
such calculations were performed in the frame of light- 
front relativistic quantum mechanics. In the present pa- 
per we calculate the " soft" part of the charge form factor 
of pion in the frame of IF RHD. We obtain a good de- 
scription of the behaviour of the form factor in the wide 
region of momentum transfers where the experimental 
data exist < Q 2 < 8(GeV/c) 2 ft}. 

We pay a special attention to the role of relativistic 
properties of quark spins for the pion structure. The 
relativistic spin rotation effect (Wigner rotation) caused 
by the summation of quark spins gives large contribution 
to the pion form factor ( 10 - 20% depending on the 
value of momentum transfer) [E4j. It is interesting that 
spin rotation effect vanishes as Q 1 — > but gives large 
contribution (~ 30%) to the pion charge radius p7| , 
which is defined by the slope of the form factor at Q 2 = 0. 

In our version of IF RHD the meson form factor asymp- 
totics at large momentum transfer is the same as in per- 
turbative QCD. The asymptotics is now determined by 
relativistic kinematics only, specifically by the relativis- 
tic effect of spin rotation, and does not depend on the 
choice of the quark wave function, that is of the quark 
interaction model. 

Our version of IF RHD approach gives, in agreement 
with the experimental data, an adequate description of 
the pion electromagnetic structure in large region of mo- 
mentum transfers. 

The pion is spinless, so the electromagnetic current ma- 
trix element has the form ( p6| ) with p c — > p n , F C (Q 2 ) — > 
F- K (Q 2 ). In the frame of composite quark model pion is 
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considered as the bound state of u— and d- quarks. We 
suppose that quark masses are equal: m u — md — M. 

To calculate in MIA the composite-system form factor 
one needs to construct first the free two-particle form 
factor ([361), ([i"o|), (|55|). Contrary to the simple model of 
the previous Section now we consider the system of two 
charged particles with spins 1/2. This gives the follow- 
ing complications. First, the Eq. ( |32| ) for the current 
operator of the free system is now transformed to the 
form: 



i£ 0) (0) = ji„ ® la 9 j 2/1 <8> A 



(82) 



Here j (1,2)11 ~ the electromagnetic currents of particles, 
^(1,2) ~ the unity operators in the one-particle state 



Hilbert spaces. The Eq.(82) can be rewritten in terms 
of matrix elements: 



{Pi ,mi;p 2 ,m 2 1 # 5 (0) \p[ , m' x ;p' 2 ,m' 2 ) = 



(p 2 ,m 2 \p' 2 , m' 2 )(pi, miljinlPi, m' x ) 



(1 



2). (83) 



Second, the matrix element of one-particle current 
contains now, contrary to ([53]) , the magnetic form factors 
of quarks as well as the charge ones. Now the parame- 
terization (the elementary-particle one following 0) is 
of the form: 

(p, m\f(0)\p\ m') 

= Wpm7 Mu p'm' Fi(Q 2 ) - UpmCT^qu Uf> m > F 2 (Q 2 ) , 



(84) 



Upyn - the Dirac bispinor, 7^ - Dirac matrix, 



Using multipole parameterization we can write the one- 
particle current matrix element in terms of Sachs form 
factors: 

G £ (Q 2 ) = F 1 (Q 2 ) + |^F 2 (g 2 ), 



G m (Q 2 )=F 1 (Q 2 ) + kF 2 (Q 2 ) , 

Fi{Q 2 )=eFi{Q 2 ), F 2 (t) = ^F 2 (Q 2 ) . (85) 

Here Ge,m - Sachs electric and magnetic form factors, 
respectively, e is the particle charge, k is the anomalous 
magnetic moment. 

It is convenient to use the canonical parameterization 
of matrix elements 132] : 



= £ (m\D*(p, P ')\m")(m"\ /i(Q 2 )^+i f 2 (Q 2 )R^\m') , 

m" 

K'^ = {p + p')^ R^e^ Xp p»p lX YP{p') . (86) 
T(p) is 4-vector of spin: 

f(p) = Mj+^M-, r (p) = (pj). 

Po + M 

The form factors fi(Q 2 ) f 2 (Q 2 ) are the electric and mag- 
netic form factors of particles. They are connected with 
Sachs form factors 153] : 



A(<3 2 ) 
h{Q 2 ) = 



2M 



y/4M 2 + Q 2 
4 

My/ 'UP + Q 2 



G E {Q 2 ) , 

G M (Q 2 ) • (87) 



Third, now the CG coefficients are of more complicated 
form. They are given by (jll]) with J = S = I = 0. Con- 
trary to the previous simple case, now the CG coefficients 
contain the Wigner rotation matrices. 

Finally, the free two-particle form factor for the system 
of two particles with spin 1/2 and quantum numbers J = 
S = I = is of the form (see also ]p0[): 



(s + af + Q 2 )Q 2 



gl q {s,Q ,s) = n c 



2y/(s-AM 2 )(s / -AM 2 ) 



( S ,QV) 



[A( S ,-QV)] 3/2 y/l + Q 2 /4M 2 
x {(« + s' + Q 2 )(G E (Q 2 ) + G|(Q 2 )) cos( Wl + lo 2 )+ 



+ — Q 2 , s')(Gli{Q 2 ) + G d M (Q 2 )) sin( Wl + u 2 ) 



Here 



£(s, Q 2 , a') = sJss'Q 2 -M 2 \(s,-Q 2 ,s') 



n c is the number of quark colours, tJi to 2 - the Wigner 
rotation parameters: 



uii = arctan ■ 



as,Q 2 ,^ 



M 



(Vs + vV) 2 + Q 2 +Via 7 (Vi+VV) 
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a(s,s')^s,Q 2 ,s') 



M(s + s' + Q 2 )a{s, s') + Vss>(AM 2 + Q 2 ) 

(89) 



L02 = arctan 



with a{s, s') = 2M + and G% M (Q 2 ) are Sachs 

form factors for quarks. The 9 - function in (J8£j) is the 
same as in (flO"|). 

An interesting effect follows from (pl|): due to the rel- 
ativistic Wigner spin rotation effect the pion charge form 
factor contains the contribution of quark magnetic form 
factors. 

The pion charge form factor can be calculated using 
(|55|), with (|8S|) for the free two-particle form factor: 



FAQ 2 ) 



d.fsdVs l V {s)gl q { Sl Q 2 : s') V {s l ). (90) 



B. The lepton decay constant of pion 

Let us calculate now the lepton decay constant of pion 
in the frame of our approach. The interest to such a cal- 
culation is threefold. First, this constant is measured in 
experiment with great accuracy ]53j], so that it can be 
a test for the model, and give the limits for parameters 
of models. Second, it is interesting to describe the elec- 
tromagnetic form factors and the weak decay constant 
in the frame of one and the same approach: the decay 
constant indirectly, through the parameters of the model, 
defines the behavior of form factors at large values of mo- 
mentum transfers Q 2 . Third, it is interesting to estimate 
relativistic effects in the lepton decay. 

The lepton decay constant f„ is defined by the 
electroweak-currcnt matrix element [h6f: 



1 



(2tt) 3 / 2 



(91) 



- 4-momentum of meson. Let us decompose the l.h.s. 
of ( |9l"l ) in the basis ( |9|). Using the explicit form of the 
meson wave function (|2lD one can obtain for (91): 



Nr. 



CG 



cVi (0^(0)1^,^(5) 



ifir Pit , 



(2tt) 3 / 2 



(92) 



As in Section II (Eq.@) one can divide the integrand in 
( |92| ) into two parts: the covariant part (smooth ordinary 
function) and the invariant part. 



Nr. 



N 



CG 



(0\jJO)W n ,V^)=iG(s)BJs) 



(2tt) 3 / 2 



(93) 



The invariant form factor G(s) is a generalized func- 
tion. In the same way as in calculating ( |54| ) of the previ- 
ous section, we now obtain the lepton decay constant of 
pion in the form 



dy/s G(s)ip(s) = fir . 



(94) 



In general, the form factor G(s) can be calculated in 
the frame of the standard model for electroweak inter- 
actions. However, in this paper we limit ourselves by 
4-fermion interaction. We take for G(s) the form factor 
which parameterizes the decay of free two-quark system: 



1 



(2tt 



,3/2 



(95) 



The explicit form (|95|) is written by analogy to (|36|) not 
taking into account the current conservation law. The 
form (|9^) is quite similar to (^lj) but instead of the con- 
stant /,r the form factor depending on invariant variables 
is written. To calculate Gq(s) let us decompose ( |95| ) in 
the one-particle basis M). Now we obtain for (Ba): 



iG (s)P^ 



(2tt) 3 / 2 



E 



mi ,m,2 ,2 C 



dpi dp 2 . , .(o) 1 -, , 

2pio 2p 2 o 



x (fii ,mi ;P2 ,m 2 \P ,\/s). 



(96) 



i c = 1, 2, 3, the sum over i c is the sum over the colours. 
The CG coefficients are known (Jl9|) . The current matrix 
element in the basis (^) can be written in the standard 
way in terms of the lepton decay current matrix element 



(2tt); 



(0|4 0) Ipi ,mi;p2,m 2 ) 
■ v(p 2 , 7713)7^(1 + 7 5 )"(Pi i m i) ■ 



(97) 



We integrate in (|96| ) in the coordinate frame with P = 0. 
Finally, we obtain: 



Go (a) = 



2V2ttP 



(po + M) 



1 - 



(pa + M) 2 



(98) 



Po 



= Vk 2 + M 2 



Substituting ( p8[ ) in the Eq.(J94|) we obtain the result 
which has the following form if written in invariant vari- 
ables: 
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f 2Mnc f a r 1 ( \ 



(99) 



Let us notice that the Eg. (|99[) coincides with that ob- 
tained in the frame of light -front dynamics How- 
ever, although all forms of RHD are unitary equivalent 
fiof , nevertheless after the physical approximations are 
made in more complicated cases the results, e.g. for 
form factors, can be different. This is possibly due to 
the fact that the unitary operators connecting different 
forms of RHD are interaction dependent JnJ and so the 
RHD forms realize one and the same approximation in 
different ways. 

Let us remark that the nonrclativistic limit of the 
Eq.(|99"|) gives the standard form in terms of coordinate 
space wave function at zero value. 



C. The results of calculations 

To calculate the electroweak structure of pion using 
©, ©, ©, (H) the following meson wave functions 
were utilized: 

1. A gaussian or harmonic oscillator (HO) wave func- 
tion 



u(k) = N HO cxp(-k 2 /2b 2 ) . 
2. A power-law (PL) wave function 

u(k) = N PL (k 2 /b 2 + l)" n , n = 2 , 3 



(100) 



(101) 



3. The wave function with linear confinement from 
Ref. [EH: 



u(r) = N T exp(-ar 3/2 - fir) , a — -vM a , 



>-T> 



(102) 



a, b - parameters of linear and Coulomb parts of potential 
respectively. 

In the Ref. ]2(J in the calculation of pion electromag- 
netic structure we supposed the quarks to be point-like. 
The results of J|(| can be considered as preliminary re- 
sults. However, one has to take into account the structure 
of constituent quarks fl55| ], in particular, the anomalous 
magnetic moment. As anomalous magnetic moments are 
connected with finite size of quark, one has to take into 
account the explicit form of quark form factors entering 
( |S8| ) and the pion charge form factor (|9(]). As in |l9| let 
us use the following forms for quark form factors: 

G q E (Q 2 ) = e q f(Q 2 ) , 



G q M {Q 2 ) = {e q + K q )f{Q 2 ) 



(103) 



the quark charge, 



the quark anoma- 



Hcrc a 

lous magnetic moment (in natural units). To obtain 
the explicit form of the function f(Q 2 ) let us consider 
the asymptotics of pion charge form factor as Q 2 — > 
oo , M ^ 0. 

To obtain the asymptotic behavior let us first make 
the asymptotic estimation of the integrals in (|h]) in the 
po int- like quark approximation (f(Q 2 ) = 1 , k = 
in^lOS) ). Omitting the details of calculation (given in 
[56[]) we write the final result for the asymptotics in the 
form: 



FAQ 2 ) 



(104) 



The asymptotics does not depend on the actual form of 
the wave function and coincides with that obtained in 



QCD. The actual form we obtain, e.g. for (100) is 



F n (Q 2 ) ~ 32^2 



[r(f)] 2 & 2 
Q 2 



(105) 



It is worth to compare the form (105) with the detailed 
QCD result Ejl: 



F n (Q 2 ) 



(106) 



If a s /7r - 0.1 then (h05|) and (|106|) coincide at b ~ 
0.1. So the asymptotics (104) is quite realistic. 

In the case of non-point-like quarks we obtain another 
asymptotics because the form factor depends upon the 
momentum transfer. It is known that QCD gives log- 
arithmic corrections to (106). To agree with this QCD 
corrected asymptotics we can, for example, choose the 
following form for f(Q 2 ): 



f q (Q 2 



l 



1 + ln(l 



(r»)QV6) 



(107) 



Here (r 2 ) is the MSR of the constituent quark which can 
be considered as the model parameter. Let us fix it (as 



in (19)) to be: {r 2 q ) ~ 0.3/M 2 . 



For the constituent quark mass in pion we use the value 
which is usually used in the calculations in RHD: M — 
0.25 GeV. 

The quark anomalous magnetic moments can be taken 
from H|: K u = 0.029 , K d = - 0.059. 



W e choose the parameters b in (100), (101) and a in 
( |102| ) in such a way as to fit the pion MSR: (r 2 ) — 
(0.432 ± 0.016 ) Fm 2 ||. We choose this way to fix 
the model parameters because the pion MSR is defined 
by the form factor at small values of Q 2 , that is the range 
where potential models work well. 
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The fit of the pion MSR gives the follo wing parameters 
of the wave functions: in the model fllO(i|) b = 0.2784 GeV, 
model JlQ]] ) at n = 2 b = 0.3394 GeV, model ( |l0l|) at n 
= 36 = 0.5150 GeV, model (|l0|) b = (4/3) a a , a s = 0.59 
at light meson mass scale, a = 0.0567 GeV 2 . 

The results of calculation are presented on Figs. 3 and 



F„ 2 (Q 2 



0.8 



0.6 




0.00 



0.04 



0.08 



0.1 2 



Q (Gev/c) 

FIG. 3. The square of the pion form factor at small values 
of momentum transfers for different models. 



The square of the pion form factor at small values of 
momentum transfers for different models (10C) - (102) is 
prese nt ed o n Fig. 3. Result s of calculation in the models 
( p00| ), ( pMl ) at n = 3 and (jl02j) coincide very closely. 

The calculations of product Q 2 F 7r ( Q 2 ) at high mo- 
mentum transfers for different models ( |100| ) - (102) are 
presented on Fig. 4. Le gend is following: 1 - harmonic os- 
cillator wave function ( [L00|) , 2 - power-law w ave f unction 
(101) at n = 2, 3 - power-law wave function (101) at n — 



3, and wave function from model with linear confinement 



(102) (these curves coincide very close ly). 



All the models for the interaction (|10C| ), (p.01|) , (|102|) 
rive a good description of the existing experimental data 



The dependence of the results on the actual model 
is much less pronounced that in the case of point-like 
quarks pC[ . 

The lepton decay constants calculated following Eq. 



( p9| ) with different wave functions hav e the following val- 
ues: U = 0.1210 GeV in the model fllPO) ), f„ = 0.1327 
GeV in the model ( p0l|) with n = 2, f„ = 0.1282 GeV in 
the model rflOl]) with n = 3, and f„ = 0.1290 GeV in the 
model (102). Let us emphasize that we have used no fit- 
ting parameters to calculate the lepton decay constant. 
Nevertheless, the obtained values are very close to the 



experimental value: 



0.6 - 



exp 



0.1307 ± 0.0005 GeV [63 



0.4 - 



0.2 - 




0.0 

4 

Q 2 (GeV/c) 2 

FIG. 4. Electromagnetic form factor, Q 2 F^(Q 2 ), at high 
momentum transfers. 

Now let us compare the numerical results for the pion 
form factor obtained in MIA (90) with that of the tradi- 
tional IA. Let us choose for the comparison, for example, 
the null-component of the current. 

To obtain the pion form factor in IA we proceed in 
the same way as while obtaining (|59|) of the preceding 
Section. Now, however, 



1) the decomposition (126 
element over the state setu 



of the IA matrix current 
is realized following (|S3"1), 



2) the parameterization of the one-particle matrix el- 
ement is given by (|8^), (instead of ([33])), 

3) the CG coefficient ( p9[ ) in (|5^) are for pion quantum 
numbers. 

Acting in the same way as while obtaining (|5^), and 
using the null-component of the current matrix element, 
we can write the pion form factor in IA in the following 
form: 



F n (Q 2 



M v \J2{2 M% + Q 2 ) 

4 4Af 2 + Q 2 yr 



2 /4M 2 



4 TheJLab new results |5£ 
our approach in ^] 



are discussed in connection with 



dy/s d 



' y/(s — 4M 2 )(s' — AM 2 ) 
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(s + s' + Q 2 ) 3 Q 2 1 



[A(s. 



V. CONCLUSION 



Let us summarize the results. 



x {(s + a' + Q 2 ) [g u e {Q 2 ) + G|(Q 2 )] cos(wi + u 2 ) 

+ jj^s,Q 2 ,s') [G U M (Q 2 ) + G d M (Q 2 )] sin( Wl +c 2 )} . 

(108) 



Here A'U = 139.5702±0.0004 MeV )53| is mass of pion. 

The normaliza tion condition F^O) = 1 is satished for 
the form factor (108) if the wave functions ( |22| ) satisfy 

To compare the numerical re sults given by the 
Eqs.© 



( |88j ) with that given by (108) let us calc ulate 
the pion form factor using the wave function ( 100 ) with 
the parameters of the calculations presented in Figs. 3 and 
4. The results are shown in the Fig. 5. The results ob- 
tained with the use of the parameterization (|55|), d88| ) 
differ essentially from that obtained without such pa- 
rameterization (108). The form factor calculated in our 
approach describes the existing experimental data ade- 
quately. 




0.6 
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0.2 



0.0 

4 

Q 2 (GeV/c) 2 

FIG. 5. Q 2 F(Q 2 ) for MIA (1) and for IA (2). Results 
of calculation with wave function (100). Parameters are the 
same as in Fig. 3. 

Let us emphasize once again that the form factor ob- 
tained in MIA does not depend on the choice of coor- 
dinate frame. This is an important advantage of our 
relativistic MIA. 



1. A new approach to the electromagnetic proper- 
ties of two-particle composite systems is developed. 
The approach is based on IF RHD. 

2. The main novel feature of this approach is the new 
method of construction of the matrix clement of 
the electroweak current operator. The electroweak 
current matrix element satisfies the relativistic co- 
variance conditions and in the case of the electro- 
magnetic current also the conservation law auto- 
matically. 

3. The method of the construction of the current op- 
erator matrix element consists of the extraction of 
the invariant part - the reduced matrix element on 
the Lorentz group (form factor) - and the covariant 
part defining the transformation properties of the 
current. The form factors contain all the dynamical 
information about transition. The properties of the 
system as well as the approximations used are for- 
mulated in terms of form factors, which in general 
have to be considered as generalized functions. 

4. The approach makes it possible to formulate rela- 
tivistic impulse approximation (modified impulse 
approximation - MIA) in such a way that the 
Lorentz-covariance of the current is ensured. In 
the electromagnetic case the current conservation 
law is ensured, too. 

5. The results of the calculations are unambiguous: 
they do not depend on the choice of the coordinate 
frame and on the choice of "good" components of 
the current as it takes place in the standard form 
of light-front dynamics. 

6. The formalism enables one to solve in part the 
problem of connection of RHD and QFT by com- 
parison of RHD with the dispersion approach. 
In this paper RHD is compared with a modified 
approach where dispersion-relation integrals over 
composite-particle mass are used. 

7. The effectiveness of the approach is demonstrated 
by the calculation of the electroweak structure of 
the pion. Our approach gives good results for the 
pion electromagnetic form factor in the whole range 
of momentum transfers available for experiments at 
present time. 
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